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Abstract 

In general the composition of Fourier integral operators (FIOs) 
need not be an FIO. Motivated by the problem of linearized seismic 
inversion in the presence of cusp caustics for the background sound 
speed, we consider FIOs whose canonical relations have certain two- 
sided cusp degeneracies, and show that the resulting compositions 
have wave-front relations in the union of the diagonal and an open 
umbrella, the simplest type of singular Lagrangian manifold. 

1 Introduction 

A fundamental problem concerning Fourier integral operators (FIOs) is that, 
outside of the standard transverse pTO] and clean intersection O [32] calculus, 
a composition of two FIOs is typically not another FIO. Describing the oper- 
ators resulting from the composition of completely general FIOs and placing 
them in a usable class, with a symbol calculus, Sobolev space estimates, and 
the possibility of constructing parametrices under suitable ellipticity con- 
ditions, is at this point a distant prospect. Some progress has been made 
for specific geometries arising in integral geometry |il6j and inverse problems 
P, 13 EH EH]. In all of these works, the compositions were shown to belong 
to an existing class, the pseudodijferential operators with singular symbols, 
which are not FIOS and were introduced to construct parametrices for op- 
erators of real or complex principal type [231 [181 126]. Such operators have 
wave front relation, i.e., the wave front set of the Schwartz kernel, contained 
in the union of the diagonal and another smooth canonical relation which 
intersects the diagonal cleanly. 
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In the present paper, motivated by a linearized inverse problem from 
seismology for the acoustic wave equation, we analyze FIOs having a special 
type of degeneracy, which we call a flat two-sided cusp, and show that their 
composition results in operators which are fundamentally more singular than 
those in works referred to above. 

FIOs with two-sided cusps, i.e., those for which the projections from the 
canonical relation both to the left and the right are cusps, arise naturally 
when considering generalized Radon transforms. They have previously been 
studied in terms of their Sobolev regularity properties; it was shown in 
[H [13] that there is a loss of 1/4 derivative compared with the nondegen- 
erate case of local canonical graphs. — > L'^ estimates for FIOs with 
one-sided cusps have also been obtained, with T*T making an appearance 
through Strichartz estimate type arguments |T2] ; however, such compositions 
have not been studied in their own right. We show that for the subclass of 
flat two-sided cusps which we formulate, the normal operator F*F is not a 
standard pseudodifferential operator with singular symbol, as it is in the case 
of fold caustics [2S1 El E]- This includes the canonical relations underlying 
the linearized forward scattering operators F in the seismology problem in 
the presence of caustics of cusp type. While there is still a large part of its 
wave front relation, WF{Kf*f)' , which is contained in the diagonal. A, the 
remaining portion, C, is now no longer a smooth canonical relation. Rather, 
C has the structure of the simplest kind of singular Lagrangian manifold, an 
open umbrella [H]. The diagonal and this open umbrella intersect in codi- 
mension one, and we quantize the variety A U C by associating to it a class 
of generalized Fourier integral operators, /"'(A U C). This is in the spirit of 
the paired Lagrangian (or J^'') distributions associated to a pair of cleaning 
intersecting smooth Lagrangians; however, rather than nondegenerate phase 
functions and product type symbols as in [211 HSl ES] , we use a combination 
of a degenerate phase function (whose gradient exhibits normal crossings) 
and standard symbols. We remark that operators with wave front relation 
in the union of A and a canonical relation having a conical singularity were 
considered by Melrose and Uhlmann [25] . 

As applied to the seismic imaging problem, the results here have nega- 
tive implications (as do the results of E] in the case of fold caustics): 
even microlocally away from A fl C, the non-pseudodifferential part of the 
normal operator has the same order as the pseudodifferential part, resulting 
in a strong, nonremovable artifact. However, our work provides a precise 
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description of the imaging artifacts resulting from cusp caustics, observed by 
Nolan and Symes [22] in a model 2D case. It would be of interest to obtain 
Sobolev estimates for operators associated with A U C, similar to the results 
in [8] for certain P''- classes, including those arising from fold caustics. 

We begin in Sec. [2] with a simple example from harmonic analysis of 
a FIO having two-sided cusps of the type that is relevant for the seismol- 
ogy problem, recall some basic singularity theory, and examine the corre- 
sponding normal operator. In Sec. [3l we describe the linearized inverse 
problem from seismology, and show that the canonical relation underly- 
ing the linearized forward scattering operator F is (i) associated with a 
two-sided cusp, but with the additional (highly nongeneric) features that 
(ii) the cusp points for the left and right projections are equal, and (iii) the im- 
ages of the cusp points on both the left and right are involutive (coisotropic). 
In Sec. m we define a general class of canonical relations having this struc- 
ture, and then derive a weak normal form, very close to the model examined 
in Sec. |2l for a general flat two-sided cusp. Finally, in Sec. [5l we use the 
oscillatory representations found in Sec. IHto analyze the composition B*A 
for two FIOs associated with such a canonical relation. We point out that 
the technique of deriving weak normal forms for degenerate FIOs and using 
them to prove composition theorems has its origins in [16] and was continued 
in [7]; different notions of weak normal forms were used to obtain estimates 

in [mis]. 

2 Example of an FIO with a flat two-sided 
cusp 

2.1 A generalized Radon transform 

We start with a simple example of a generalized Radon transform that in- 
corporates the features of the linearized seismic inversion problem in which 
we are ultimately interested. In fact, we will later show that perturbations 
of this example serve as weak normal forms for the general class of operators 
with fiat two-sided cusps. 

The operator averaging over translates in of a curve 7(t) such that 

7, 7, 7 , 7''''^ are linearly independent (1) 
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is an FIO having a canonical relation which is a two-sided cusp [12] , but the 
same holds for the curve 7(t) = (t, t^, t^) in [131 IE] , and this is the model 
we will use. 

Thus, consider the generalized Radon transform Rq : V{M.^) — > V{M.^), 



where x ^ C^(]R) is a fixed cutoff function. The associated canonical re- 
lation, the conormal bundle Cq = N*Z' C {T*R^ \ O) x {T*R^ \ O), where 
Z = supp (Krq) C M^xM^ and (x, ^; y, r])' = (x, ^; y, —rj) is standard notation 
for the twist map, is 



We will show that Co has several properties, which remarkably also hold in 
the totally unrelated seismic imaging problem: (i) the projections both to the 
left and right, hl, hr : Cq — > T*M^\0, have (Whitney) cusp degeneracies; (ii) 
the cusp points for the two projections are the same. Si i(7ri) = Si i(7r^) := 
Si i; and (iii) the images of the cusp points, 7ri(Si i) and 7r/j(Si i), are 
coisotropic (involutive) submanifolds of T*]R^. Conditions (ii) and (iii) are 
unstable and quite special among two-sided cusps, i.e., canonical relations 
satisfying (i); see the Remarks in Sec. 12.51 below. 

2.2 Singularity classes 

We first recall some basic facts about cusps and refer to [531 EZl [ID] for 
more details. Let / : — be a smooth function. We say that / 
drops rank simply at p if rank {df)p = N — 1 and if {d{det df))^ ^ 0, so 
that Si(/) := {x G M.^ : det{df{x)) = 0}, the corank one critical set of /, 
is a smooth hypersurface near p. If kerd/p ^ TpSi(/), then / has a fold 
singularity, and is an immersion. Considering the more degenerate 

case when kei dfp C TpSi(/), one may choose a nonzero vector field v along 
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xi, X2, X3, -2(xi - yi)92 - 4(xi - yi)^93, 02, 03; 
yi, y2, Vs, -2(xi - yi)92 - 4(xi - yi)^93, O2, : 
X2-y2- {xi - yif = x3-y3- (xi - yi)^ = 



(3) 
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Si(/) such that v G kerrf/p, so that v is tangent to at p. Let g he a. 

smooth function such that 5'|si(/) = and dgp ^ 0, e.g., g = det df. Thus, 
dg{v) has a zero at p. 

Definition 2.1. / has a (Whitney) cusp at p if dg{v) has a simple zero at p. 

the c«sj> sei of /, is then codimension 2, and /Isj is an im- 
mersion. One may use adapted coordinates to clarify this. These are lo- 
cal coordinates such that /(xi, X2, . . . , xtv) = {xi,X2, ■ ■ ■ ,XN-i,h{x)) and 
h{0) = [27]. Then, = {x : ^(0) = 0}, and / has a cusp singu- 

larity at iff ||(0) = 0, ||(0) ^0 and rank = 2. 

These conditions are adapted coordinate- independent, and the notion of a 
cusp makes sense for any smooth mapping between A^- dimensional man- 
ifolds. In suitable smooth coordinates on the domain and range spaces, 
any map with a cusp singularity can be put into the local normal form, 

/(Xi, X2, . . . Xn) = {Xi, X2, . . . , Xn-1, Xn-iXn + 

For the canonical relation Cq, the projection to the left, txl '■ Cq — > T*M.^, 
is nL{xi,X2,X3,yi,92,93) = {xi,X2,X3, -2(xi - yi)92 - 4(xi - yif 63, 92,93); 
hence, letting a = 2^^2 + 12(xi — yi)^6'3, f3 = — 2(xi— yi) and 7 = — 4(xi— 
one has 

" 1 " 
1 
1 
—a a p 7 
1 
1 

which has det diiL = 292 + 12(xi — yi)^93 and kervrj^ = M 
{2^2 + 12(xi - y^y93 = 0}. Since ^(2^2 + 12(xi - yi)^) 
and ^(2^2 + 12(xi - 2/1)^^3) = 24^3 0, the cusp set is 

Si,i(7r,.) = {2^^2 + I2{x, - yi)% = xi - y^ = 0} = {92 = x, ~ y, = 0}. 

Noting that 26^2 + 12(xi — yiY93 and Xi — yi have linearly independent gra- 
dients, we see that ttl has a cusp singularity. 
Similarly, the projection to the right is 

TTnix, yi, 92, 93)={yi, X2-(a;i-?/i)^ X3-{xi-yi)'^, -2{xi-yi)92-4:{xi-yi)^93, 6*2, 93), 



■ ^ at Si(7rL) = 
= -24(xi-yi)^^3 
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so that detdnji = —26*2 — 12(xi — 7/1)^^3, andkervr/j = at Ei(7rij) = {26^2 + 
12{xi-yi)% = 0}. One also has £-^{292 + 12{xi-yi)^97) = 24{xi-yi)93 and 
^(2^2 + 12(xi - yife^) = 24^3 ^ 0, so that Si,i(7r^) = {6^ = xi-y, = 0}, 
and thus ttji also has a cusp singularity. 

Note that T,i{nji) = ^1(7?^) and T^i^i^nn) = ^1^1(7?^) =: Si^i. The first 
of these is true for any canonical relation [19\, but the second is a very 
strong condition; indeed, for general two-sided cusps, there is no relationship 
between ^11(71^^) and Ei i(7rj:j). Furthermore, for the images of these sets, 
one has 

vrL(Si) = {^2^ = -f vr«(Si) = {4 = -^vhs}, 

with 7ri(Si_i) = {^1 = 6 = 0}, 7rR(Si,i) = {771 = 772 = 0}, resp., their 
cuspidal edges. 7ri(Si 1) and 7r^(Si 1) are codimension two coisotropic (or 
involutive) submanifolds of T*]R^\0 [19j; again, this is a nongeneric situation. 

Consider next the composition of canonical relations, 
Co°Co = {{x,^;y,r]) : 3(z,C) s.t. {x,^]z,C) e and {z,C,y,v) e Co}. 
For such a triple {x,C,; z, y, rj), one obtains the following equations: 

Z2-X2 = {Zi - Xif, Z3-X3 = {zi - Xif , 

^1 = = -2(2:1 - Xi)C2 - 4:{zi - xifCs, 6 = C2, 6 = Cs, 
Z2 - y2 = (zi - yi)"^ , z^-y^ = {zi-yiY, 

Vi = Ci = -2(2:1 - yi)C2 - 4(zi - yifCs, V2 = C2, V3 = Ca- 
using 2(^1 - xi)^2 + 4(2:1 - xi)3^3 = 2(2:1 - yi)C,2 + 4(2:1 - yifC.3, after sim- 
plification one obtains: 

{yi - a;i)6 + 2[(zi - x,f - {z^ - yi)% = 0, 

(2/1 - xi)[^2 + 2{3zl - 3zi{xi + yi) + xl + yl + xiyi)^:^] = 



6 



It follows that the contribution to Cq o Cq from {yi — Xi = 0} is contained 
in A, and that from {yi — Xi ^ 0} is contained in Cq, where 

Cq = ^^{xi,X2,X3,-2{zi-Xi)92-4:{zi-Xif 93,92,93; 

yi, y2, Vs, -2(2:1 - xi)92 - A{zi - 0:1)^^3, ^2, ^3) : (4) 
X e M^ vi, zi G M, ^3 e M\0, y2 = X2 + {yi - xi){2z^ -x^- yi), 
y3 = X3 + {yi - xi){2zi - xi - yi){{zi - xif + {zi - yif), 

6*2 = -2(3^1 - 3^1 (xi + yi) +xl + yl + Xiyi)93 | 

Notice that ii = Tji = 46'3(zi — xi){zi — yi)(2zi — xi — yi) and that Co 
intersects A in codimension one, at {xi — yi = 0}. 
The parametrization of Cq in (jlj) is a map 

■■ ^'„..,.3,.i,.i X (^^3 \ 0) (r*M3 \ 0) X (T*M3 \ 0). 

One easily sees that T is singular at S := Si(T) = {xi — yi = xi — zi = Qi} , 
where kerciT = M ■ which is ^ TSi. Thus, as discussed below, Cq is 
an open umbrella, exhibiting the simplest kind of singularity of a Lagrangian 
manifold. 



2.3 Open umbrellas 

For maps between manifolds of the same dimension, a fold is a singularity 
of the type Si^ and a cusp is a singularity of type 5'i^i^o [Ml EZl ID]- From 
larger dimensional spaces to smaller ones, Si^q maps are submersions with 
folds, while in the opposite direction they are referred to variously as cross 
caps or Whitney- Cay ley umbrellas [T0l[9]. In the lowest possible dimensions, 
if g : ^ has a cross cap singularity then, in suitable local coordi- 
nates, {u, V, w) = g{x, y) = (x^, y, xy), and its image is the algebraic surface 
{w^ = Mv^}, the Whitney-Cayley umbrella. This is the simplest type of non- 
immersed surface singularity in three dimensions. It is an immersion away 
from the origin; it is actually an embedding on {y 7^ 0}; and along {y = 0}, 
it is folded and hence 2-1. 

Now, by adding one dimension to the range space, one may simulta- 
neously both unfold the closed umbrella, making the parametrization 1-1 
away from the origin, and make it a Lagrangian manifold with singularity in 
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M^. The open (or unfolded) umbrella is the map U : ^ R^, U{x,y) = 
(x^, y, xy, (and its image). We have U*u = 0, where u is the symplectic 
form on ~ T*M^. Hence, the unfolded umbrella is a Lagrangian inclusion, 
i.e., a variety which is a smoothly immersed Lagrangian manifold away from 
its singular points [9]. To find a similar conic model, one needs to look in 
T*i?^\0; for an unknown function /(t, s), consider 

= { t, /(t, s); ste, -t^e, e) e T*R^ \ O : ^ ^ O }. 

3 

Then, Ai is a Lagrangian if / is chosen so that A dx\A-^ = 0. Since 
d^Adx = d{st9)Ad{t'^) + d{lt^9)Ads + deAdf = dO A{2st'^dt + lt^ds) + de Adf , 
this holds if df = —2st'^dt — ^t^ds = —d{^t^s); using / = —^t^s yields the 
conic Lagrangian, 

2 2 

Ai = { s, — t^s; ste, -t^e, e) : t,s eR,e eR\o}, 

3 3 

exhibiting an open umbrella singularity along {(0, 0, 0; 0, 0, ^3) : ^3 7^ 0}. 

Observe that Ai may be parametrized by a degenerate phase function. 
Letting (p{x, Oi, 6*2, r, rf) = [xi - {^f )Oi + [xs + |(p^a;2)6'2, one has 

T 2 T 

doj^if = xi - {-f , d02(f = X3 + -{-Yx2 

T] 6 7] 

drv = -(-2^ + 21^x2), d,^ = -C-n-'i- + 2-^x2), 

7] r] 7] rj 7] 7] 7] 7] 

SO that (p is degenerate, with both dr^p and djj^p having normal crossings: 
dr^p = iff r = or —2^ + 2|^X2 = 0, and d^y? = on same two sets. 
For the second case, solving for 61 = ^X2d2 we see that ip parametrizes Ai, 
while for the first case, (p also parametrizes Aq := A^*{xi = X3 = 0}. Ag 
and Ai intersect cleanly in codimension one, except at the singular set of 
Ai, which is contained in the intersection. ip> simultaneously parametrizes 
Aq U Ai, although it is not a multiphase function in the sense of [26], used to 
parametrize a pair of cleanly intersecting smooth Lagrangians. The class of 
generalized FIOs we define in Sec. 15.11 has similar features. 

For the general notion of an open umbrella, let (M, uj) be a smooth sym- 
plectic manifold of dimension 2n. 
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Definition 2.2. S" C M is an open umbrella if it is the image of a map 
ip : -R" — )■ M such that dijj drops rank by one simply at Si, of codim two; 
ker TT,i ; and Range (dipp) is Lagrangian for p e i?"\E. 

The class of open umbrellas is structurally stable under smooth pertur- 
bations in the class of Lagrangian inclusions [H]. As checked above, the 
parametrization T in (jl]) satisfies the conditions in the definition, so that Cq 
is an open umbrella. 

2.4 Composition for the model fiat two sided cusp 

We now want to show that the composition A* A for Fourier integral op- 
erators, such as Rq, associated with the model flat two-sided cusp Cq in 
([3]) results in operators with wave front relation in the union of the diago- 
nal A and the open umbrella Co from dl). Let A e J™(]R^,]R^; Co), so that 
Af{x) = J e'^'^-'y''^a{x,y,e)fiy)dedy with (l){x,y, 6) = {x2-y2-{xi-yif))e2 

+ (3^3 ~ Us ~ i^i ~ 1/1)^)^3 and a G S^^^, the standard Hormander class of 
symbols. We consider the normal operator formed by the composition, 

A*Af{x) = j e'^'^^''y^'^^-^^'^''^^^^aiz,y,ri)a{z,x,9)fiy)dydzd9dri. 

We have 

<l>{z,y,v) -Hz,x,6) = {z2 - y2 - {zi - yiy))v2 + {z3 - y3 - {zi - yi)^)m 

-{Z2 -X2- {Zi - Xif))92 - (^3 -X3- {Zi - Xi)'^)93. 

After a stationary phase in {z2,'r]2) and (2:3, r/3) in the integral for Kyi*A, we 
obtain: 

(l>iz, y, v)-(l>iz, X, 9) = ix2-y2+{zi-xif-{zi-yif)92+{xs-y3+izi-xi)'^-{zi-yi)'^)93 

= {x2-y2)92+{x3-y3)03+{yi-xi)[i2zi-xi-yi)92+i2zi-yi-xi){{zi-XiY+{zi-yif)93] 
Let zi = ^, and then consider the change of variable, 

ei = -(2^ -X,- y,)92 - (2^ -y,- x^){d - x^)' + d- y,f)93, 
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for which 

When ^ = then 2^ = xi + yi, = and 02 = -^^^^^3- Let 
^(^1, ^2) = (ii^2 + 1^1, ^2) be a normal form for a two-dimensional map with 
a cusp singularity. To understand the push forward T^u of a distribution u 
under the map T, note that the kernel of is 

KT^{z,t) = 6{z-T{t)) = J e'[(^i-*i*2-^*?)^i+(^2-*2)^2]l((T)d(T, 

so that 

Returning now to A* A, letting zi — ^ and Z2 — ^ + — yi)^, the 
phase function of A* A becomes 

y, 0, a) = {xi - yi)9i + {x2 - ^2)6*2 + {^3 - ^3)613 

+(| - tit2 - ltl)a^ + (| + - y.f - ^2)^2 

and the amplitude becomes: a x Next, we perform stationary phase, 

first in t2, 02 and then in d\^a\. We have 

~ ~ \ 

^3 ^ 

The Hessian is 1, the phase function is = (a; — y)9 + — ~ + 
l{xi - yi)^])(7i, and 

~ (Ti ,T^1 Iq /^2 1/ \2\ 

dei0 = - yi + ^, = — - -i^ - + -(xi - yi) ). 

(73 (73 ^ "3 ^ 

The Hessian is ^. Finally, we introduce the variable t = tiO^ and obtain 
X2-y2 + - yi)^ 02 + (^X3 -y3 + ^{xi - yi){jf + ^^(^1 - 
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with respect to which the Schwartz kernel has the representation 
KA*A{x,y) = J e^hd92d93dT, b e S^"". 



The wave front relation satisfies 

T^2 IrV 3 
^3 



WF{KA*Ay C {(xi,X2,X3,^ + ^^ + ^(a;i-|/i)V,^2,^3; 



c?e20 = 2:2 - Z/2 + - yi) = 0, 

(^3 

deJ> = X3 - 1/3 - ^(^1 ~ ^1) ~ ^(^1 ~ ^1) = 0' 



(ir0 = ^(a^^i - vi) + ^^(^1 - 2/1) + ^(^1 - y^f = o}- 

From the last relation, we see that the phase function is degenerate, with the 
critical set in the phase variables having a normal crossing, being the union 
of two transverse surfaces, corresponding to Xi — t/i = or ^ + |^ + \ {xi — 

Uif = 0. The points where Xi = yi give rise to the diagonal A, but with the 
parametrization actually being a cusp: 

{x,T,e) — > (^'-^ + 2 ^'^2, ^3; ^'-^+ 2 ^'^2,^3). 

On the other hand, the points where ^ + |^ + — ViY — contribute 
to WF{A*A) the set, parametrized by x e R^, yi, r e R, ^3 7^ 0, 

^0 = { {x, y, v) ■ y2 = X2 + ^{xi - yi), ^3 = 773 = ^3, (5) 

Vs = ^3 + 2^^'^^"^^'^ +2^^^i~^i)' 

6 = ?7i = -;72 + i^i - yi)\ 

^3 



6 = V2 ^ -^{xi - yi)% - 
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which is the image of p(xi, Xa, Xg, yi, 63, r) = (x, ^; y, r]) with y2, y^, r]i,V2, V3, 
$,1,^2, ^3 given by the above relations, p satisfies the conditions of Def. 12.21 dp 
drops rank simply by one at S = {xi—yi = r = 0} and ker dp = M-^ ^ TS. 

Hence Cq is an open umbrella. 

Next, we compute the principal symbols on A and Co away from their 
intersection. We have Crit^ = CritA U Crit^^ = {{x,y,T,92,93) \dr4> = 

0,dg^(f) = 0,dg^(f) = 0}. Using Hormander's formula [IH], cr = a{E^)^ where 
= |— D(iJf~l~ local coordinates on Crit^. On Critj\, local 

coordinates are (x, r, 6*9, 6*3) and Ej = I — ' ' ' 'ar'aga'aga i-i _ _^ 

D{x y T ^ 

while, on Critp^ , local coordinates are (x, Vi, r, 63) and E2 = I-fu- — ' ' fl-^ ' ' 903 1 

' Co' ^ ' '^^ <P I D{x,yi,T,02,93,y2,y3,S2) ' 

i ^^e^^ )~^- Thus, both principal symbols are singular as one approaches 

A n Co, behaving as 5~^, where 5 is the distance on either A or Cq to 
the intersection. 



2.5 Flat two-sided cusps 

Motivated by the previous example, we now define the class of canonical 
relations for which we will establish a composition calculus. 

Definition 2.3. If X and Y are manifolds of dimension n > 3, then a 
canonical relation C C {T*X \ 0) x [T*Y \ 0) is a flat two-sided cusp if 

(i) both ttl : C — )■ T*X and ttj? : C — > T*Y have at most cusp singularities; 

(ii) the left- and right-cusp points are equal: Si i(7ri) = ^11(7?/?) =: Ei i; 
and 

(iii) 7ri(Si 1) C T*X and 7r^(Si 1) C T*Y are coisotropic (involutive) and 
nonradial. 

Here, as usual, at most cusp means either cusps or folds, which are un- 
avoidable in the neighborhood of a cusp. Nonradial, meaning that the the 
radial vector field Yl does not lie in the two-dimensional Hamiltonian 
foliation of the codimension two coisotropic submanifold, is a standard tech- 
nical assumption needed to apply the homogeneous Darboux theorem in 
Sec. H 

Other examples of fiat two-sided cusps coming from generalized Radon 
transforms, verified as with the example above, are those for translates of 
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curves in satisfying ([T]), and for the translates in of Menn's surface pQ, 
Ci = N*{xs -y3 = (xi - yif{x2 - ^2) - {.X2 - 1/2)^}', 

or any other surface in whose Gauss map has a cusp. Computation also 
shows that the composition C\o Ci = A U Ci, with Ci similar to Co as in 
([5]) above. 

Remarks. 1. If (i) and (ii) in Def. 12.31 are satisfied, the structure of 
|^T*x|T7rL(Si,i) is the same as Ci;T*y |T7r_R(Si,i), since both are equivalent to the 
restriction to Si 1 C C of the degenerate symplectic form ujc '■= t^l^jJt*x = 
vr^Wj-y. Thus, 7ri(Si 1) is involutive iff 7r^(Si 1) is, and one can see that this 
holds iff the two-dimensional ker(a;(7|si) is simply tangent to Si at Si 1. This 
is in contrast to the situation when the image of the cusp set is symplectic 
(strongly noninvolutive), in which case uc is a folded symplectic form and 
ker(a;(7|si) ^ ^i everywhere [22] . 

2. Condition (i) is stable under small perturbations, since cusps are struc- 
turally stable [in]. However, given (i), condition (ii) is unstable and atypical. 
To see this concretely, consider (possibly) nontranslationally invariant fam- 
ilies of curves {'^x}x&s^^ using the framework of ^ as analyzed in terms of 
FIOs with cusps in For vector fields X, F, Z, W on M^, with X ^ 0, let 

7^(t) = exp^itX + t^Y + t^Z + t'^W). 

The associated generalized Radon transform is in J~5(]R^, M^; C*), with C a 
canonical relation for which, by [121 Prop. 6.1], vr^ (resp., tir) is at most a 
cusp if 

X,Y, Z ±-[X,Y],W ± -[X, Z] + ^[X, [X, Y]] are linearly independent. 
6 4 24 

From the calculations in [12], one sees that a necessary condition for Si^i(7r2,) = 
Si,i(7ri?)isthat (Z+i[X,F])-(Z-|[X,F]) G span {X, y}, i.e., span {X, y } 
is an integrable distribution of 2-planes, which is generically not the case. 

Nevertheless, as seen in the next section, flat two-sided cusps arise natu- 
rally in a quite different setting without translation invariance or integrability. 
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3 Linearized seismic inversion with cusp caus- 
tics 



We consider linearized seismic inversion for the single source data set, under 
the assumption that the background sound speed has caustics of at most 
cusp type, i.e., either folds or cusps. We show that the linearized forward 
scattering operator F is a Fourier integral operator associated with a flat 
two-sided cusp, C . We only briefly describe this problem, referring to [291 
|30| [28| ini [7] for more motivation and details. 

Acoustic waves are generated at the surface of the earth, scatter off het- 
erogeneities in the subsurface and return to the surface, where measurements 
of the pressure field are used to reconstruct an image of the subsurface. The 
model for the scattered waves is given by the acoustic wave equation, 

1 d'^ 

-^j^{x,t)-^p{x^) = Kx-s)6{t), p(x,t) = 0, t<0, (6) 

where x E Y = {x E R^,Xs > 0}, representing the Earth, p{x,t) is the 
pressure field resulting from a pulse at the source s at time t = 0, and c{x) 
is the unknown sound speed field. 

To make this nonlinear inverse problem tractable, one considers a lin- 
earized operator F which maps singular perturbations of a smooth back- 
ground sound speed in the subsurface, assumed known, to perturbations of 
the resulting pressure field at the surface. Thus, the linearization consists in 
assuming c to be of the form c = Co + 6c and the resulting p = Po + Sp, where 
Co is a smooth known background field. The formal linearization of ([6]) is 

1 d'^Sp 2 d'^po 

nco(M := ^rT^Tr(^'^) - ^"^pI^'^) = — ■ ^c(x), 6p = o, t<o, 

where Pq, the Green's function for Dcq. The linearized scattering operator 
is F : 5c — )■ 5p\qyx{o,t)- Under mild technical assumptions, F is an FIO 
associated with a canonical relation C C {T*X \ 0) x {T*Y \ 0), where 
X = dY X (0,T) is the data space [3ll [201 |30]. The goal is then to left- 
invert F; standard techniques suggest studying left invertibility of the normal 
operator, = F*F. 

Let H{x,^) = |(co(a;)^^ — be the Hamiltonian associated to Cq, and 
As the image of T*]R'^ \ under the bicharacteristic flow associated to H, 
which is a Lagrangian submanifold of T*]R^ \ 0. A caustic is a singularity of 
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the spatial projection vr : — )■ F. It is known that the caustics exhibited 
for generic soundspeeds Cq are the same as generic Lagrangian singularities, 
i.e., folds, cusps, swallowtails, etc. [T7]. 

If Co is such that only fold caustics occur, it was shown by Nolan [28] 
that F G I^{X,Y] C) with C a two-sided fold. Furthermore, the Schwartz 
kernel of the operator F*F belongs to a class of distributions associated 
to two cleanly intersecting Lagrangians in {T*Y \ 0) x [T*Y \ 0), with the 
corresponding canonical relations being the diagonal A and another folding 
canonical relation [281 E]- 

The next caustics to consider are those of cusp type, meaning that the 
only singularities of the spatial projection vr : — )• F are folds or cusps. A 
cusp caustic is already present in the 2D example of [22] exhibiting loss of 
regularity for F. To analyze the geometry of the canonical relation C in the 
presence of cusp caustics, we make use of the description of A^ in [28j . It can 
be parametrized by tine, the time travelled by the incident ray, and the takeoff 
direction (^1,^2,^3) ^ S^. We can change these coordinates to {xi,X2,P3)- 
Hence on Cs, X3 = /(xi,X2,P3) and (^1,^2) = igi{xi, X2,P3), g2ixi, X2,P3)) 
and As is the graph of a function VG{xi, X2,P3) which means that ^ = 
g,^ 9G ^ &G ^ f_ rj^Y^e^ |a = |/ |22 = |1 |22 ^ In this'^new 

0x2 ' ops ops axi' ops 0x2' 0x1 0x2 

setting, 7r(xi,X2,P3) = (a^i, 2:2, /(xi, X2,P3)) det dir = and cusp caustics 



occur when 



— — = ^—15- = 0, -;—iT 0, and (Vt: — , Vtt^I are linearly independent. (8) 

Op3 Opl Op'l ^ Ops Op^ ' 

Next, we parametrize the canonical relation C of F in terms of Xi,X2 and 
P3 ; (ai,a2, a/1 — the take off direction of the reflected ray, writing 

a = («!, 02); and r, the variable dual to time. 

C = {{ri{-),r2{-),tinc{-) +tref{-),pi{-),p2{-),r; Xi,X2,f{-), 

-r(co + gi{-)), -r(co i(-)a2 + (72(-)), -r(co i( Vl " I"!' + Pa))}, 

where 

/(■) = fixi,X2,P3); rj(-) = rj(xi,X2,/(xi,X2,j93),a),J = 1,2; 

tinc{-) = tinc{Xi,X2,P3); trefi') = trefi^l, X2, f {Xi, X2-, P3) ■, Oi)'i 

Pji-) = Pi(a;i, a;2,/(a;i, 0:2,^3), a), j = 1,2; gj{-) = gj{xi,X2,P3),j = 1,2; 
and Cq\-) = Co ^(a;i,X2, /(■))• 
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We have 7rfl(xi, X2 ,P3, «i, « 2, r) = {xi,X2, /(■), -^(cq ^(■)"i+fl'i(-))> -^(^o ^(■)«2+ 
g2i-)), -T"(co^(-)\/l - l"P +P3) and thus 



■ 1 











1 







M. 


dL 


dxi 


8x2 




A 


B 


C 


D 


E 


F 


G 


H 


I 




































-Co\-)V^- 







P3 



where A 



a(-r(c-^(-)ai+gl(.))) . ^ ^ a(-r(co^(.)«i+gi(.))) , 
dxi ' 9x2 ' 

^ ^ d{-Tic-H-)<^2+92i-))) . ^ ^ 0(-r(c-l(.)«2+g2(-))) . ^ ^ 



9(-^(^0''(-)Vl-|^P+P3)). 



dxi 



H 



8x2 

d(-T{c-\-),/l-\a\2+P3)) ^ 
8x2 ' 



and / 



8(-r(co^(.)ai+gi(.))) . 
8{-r{c-\-)a2+92{-))) . Q ^ 
~r{co\-)y/l-\a\2+P3)) 



8p3 



Hence detdirn = ^ and the kernel at Si = {detdirn = 0} is spanned by 
a vector vi which is jf- plus a linear combination of \jf-, if-, f-}, because 

8p3 ^ 8ai ' oa2 ' 8t J ' 

the matrix 

-rco^ Co^(-)ai + 

-rcn^ Co^(-)a2 + fi'2(-) 



rc,, a\ 





^"2 1 
2 '^0 



1 



P3 



is nondegenerate by the calculations of [28]. From the cusp conditions 
one can see that ttr has a cusp singularity. 

Similarly, 7ri(xi, X2, ai, 0^2,^3, ^) = ('"i, ^^2, W + W' Pi' P2, 1") and 



9ri I 8ri df 

8x1 8x3 8x1 

8r2 _|_ 8r2 df 

8x1 8x3 8x1 

8tref I dti. 



8x1 



+ 



8x1 

8pi I dpi df 
8x1 8x3 8x1 
8p2 _|_ dp2 df 
8x\ 8x3 8x1 



8ri _|_ dri 8f 

8x2 8x3 8x2 

8r2 _|_ dr2 8f 

8x2 8x3 8x2 

9tref _|_ 8tinc 



8x2 
8pi 8f 



8x2 
dpi _|_ 

dx2 dx3 8x2 
8p2 _|_ 8p2 8f 
8x2 8x3 8x2 





8ri 
8a 1 
8r2 
8ai 

8tref 

8ai 
8pi 
8a 1 

8p2 

8a 1 





dri 
80.2 
dr2 
80.2 

dtref 

802 



dt, 



8ri 8f 
8x3 8p3 
8t2 df 

8x3 8p3 
ef 8f I 8t 

I': 



+ 



8x3 8p3 8p3 

8pi 8pi 8f 

80.2 8x3 8V3 

8p2 8p2 df 

80.2 8x3 8p3 











1 



From 128] we have that 



8ti, 



8p3 



Q ^''1 
' 8x3 



Q §12. 



^ 

' 8x3 



^ 



dt 



ref 



8x3 



< 
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0, and that the matrix 



dri 


dri 


dri 


dri 


dxi 


8x2 


dai 


da2 


dr2 


ar2 


dr2 


dr2 


dxi 


8x2 


dai 


da2 


dpi 


8pi 


dpi 


dpi 


dxi 


8x2 


dai 


da2 


dp2 


8p2 


dp2 


dp2 


dxi 


8x2 


dai 


da2 



is nondegenerate. Thus, detdiTL = kervTi = ^ and as before one sees 
that hl has also a cusp singularity. 

Now we will show that the images of the cusp points are involutive. By 
Remark 1 in Sec. 12.51 it is enough to examine the degenerate symplectic form 
uc = vr|j,co'T*y on C and its kernel, equal to {v : ijJc{v,w) = 0,Ww G TC}. 
We have from ([8]) that Si = {{x,^);fp^ = 0} and Si i = {{x,C,);fp^ = 
fp3P3 — 0}- Since wcIei drops rank by 2 and kerdvr^ = M ■ fi, one has 
keiuc = spanjfi, ^2}, where V2 is a vector, without ^ component, that we 
now determine. At Si, one has 



'^clsi = (co^ai +gi + (cg Vl - \a\'^ + ps) f^,){dT A dxi) 
+ {co^a2 + g2 + (cq^VI - +P3)fx2)idT Adx2) 



+ TCn 



'I - a 



-Jxi){dai A dxi 



TCn 



a2 



-ai 



\a\ 



-Jxiidoii A dx2) + rcg ^ 



-«2 



1 — a 



-Jxi{da2 A dxi 



fx2){doi2 A dx2 



+ '^{dx2Co - dxiCQ ^012 + (9^2Co - - c^xiCo Vl - \(y\^fx-^{dx2 A dxi) 

+ 2Tfx^{dp3 Adxi) + 2Tfx^{dp3 Adx2). 

We need to look at the 5x5 skew symmetric matrix of uc corresponding to 
Xi,X2,ai,a2,T: 





—a 



ai 



TC. 



1 —02 f 

J Xl 



fxi, 



yiT^ 

b 



TC, 



a 


-1 -a 



= f 



TC,\1 



y/T^^'"'^ 
a2 f 

J X2 



C 



-k -k —b 

-k -k —C 









where a is the coefficient of dxi A dx2 , b is the coefficient of dr A dxi and c 
is the coefficient of dr A dxo- f is a combination of 



d 8 



8 8 8 



8x1 ' dx2 ' 8ai ' 8a2 ' dr ' 



note 
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that 



1 - 



fxi 



fxi 



-Ql f 



X2 



1 + Va/1 



which means that fo does not contain any jp-, j^- terms, but only 73^, jf-, 
This shows that kercu^ is simply tangent to Si at Si^i, and thus 7rj:j(Ei 1) is 
involutive; by Remark 1, so is 71^(21^1). 



4 Weak normal forms for flat two-sided cusps 

Next, we show that any flat two-sided cusp can be prepared, by applica- 
tion of suitable canonical transformations on the left and right, so as to be 
parametrized by a phase function similar to that for R^'m. Recall that 
one can conjugate any two-sided fold to a single normal form [221 [23], but, as 
in [7] , we will need to work with merely approximate normal forms. For 
simplicity, we restrict ourselves to the three dimensional setting of interest 
for the seismic problem. 

Assume dim (X) = dim (y) = 3 and suppose C C T*X x T*Y is 
a flat two-sided cusp. Let G Si 1. We claim that there exist local 
canonical coordinates (a;,^) on T*X near := 7rx,(c*') and [y^r]) on 

T*Y near {y^,ri^) := T^nic^) such that, as for the model Cq in Sec. |2l 
{xi,x',yi,ri') := (xi, X2, X3, ?/i, 772, r^s) form local coordinates on C near c°. 
In fact, since and dn^f drop rank by 1 at c°, there exist (see [12]) sym- 
plectic decompositions T(^o,50)(T*X) = 14 © V, T^y0^^o){T*Y) = Wi ® W\ 
with dim {Vi) = dim {Wij = 2, dim (V) = dim {wj = 4, and 

T,oC = Gr{x) © {Al X (0)) © ((0) x A^), 

where C Vi, Aji C Wi are Lagrangian, i.e., one-dimensional, subspaces 
and X ^ Sp(y' ,W'). Thus, in suitable hnear symplectic coordinates, A^ = 
{^^ = 0}, Ar = {rii = 0} in V^i ~ VTi ~ T*R, and x = / e Sp{T*R'^,T*M.'^), 
so that {xi,yi,x',T]') are coordinates on T^oC. By a standard application of 
the homogeneous Darboux theorem [191 122] , there are local canonical coor- 
dinates {x, ^, y, rj) on T*X x T*F such that c° = (0, 63*; 0, 63*) and (x, yi, -q') 
form local coordinates on C . 

We further prepare C by noting that if vr, vl are kernel vector fields 
for ttr, Til, i.e., nonzero vector fields along Si generating ker((i7r^), ker^diiL), 
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resp., then vl is a linear combination of dy^^, dr)2, drj^. Thus, dT:R{vi) is 
a vector field along 7r^(Eij), nontangent to Si i since ttr is a cusp. Us- 
ing homogeneous Darboux and the nonradiality of 7ri^(Si^i), one can as- 
sume that vl = dy^ at Si^i, txr{^i,i) = {rji = ?72 = 0} and diiRiTC) = 
span{TnR{T,i^i) , d7iR{vL)} = {drji = 0} along Ei^i. Working similarly on the 
left, one can assume vr — d^^, T^LiJ^\,i) — {^i = ^2 = 0} and dT^iiTC) — 
{d^i = 0} along Ei i. 

Furthermore, near one can assume that (i) Si = {ri2 = f{x, yi, rj^)}, and 
(ii) Ei,i = {r]2 = f{x,yi,r)3), yi = g{x,r]3)} for functions f,g^ homogeneous 
of degrees 1,0, resp., and with dx^g ^ 0. To see (i), let f{x,yi,r]') be a 
defining function for Ei, homogeneous of degree 1; by the implicit function 
theorem it suffices to show that dr^^ f{c^) ^ 0. Euler's identity imphes that 
Tj'-dri'f = at Ei; since ?72 = 0, 773 7^ at c°, this yields c?^,^/ = 0. Since hl, ttr 
are cusps, the kernel vector fields vl, vr G TEi at Ei^i, thus at c°, and hence 
c?xi/(c°) = dy^f{c^) = 0. If dri^f = 0, the only nonzero components of df{c^) 
would be dx'f, and this contradicts the fact that dniiT'Li) = d7TL{TT,i^i) = 
TttlCZi^i) = {d^i = d^2 = 0}. Hence, d^j^fic^) 7^ 0, and thus Ei can be 
described as {rj2 — /(a:, r^s)} for some /. Note that since vr,vl £ TEi 
at Ei^i, one has 9xi/|si.i = 5i/i/|ei,i = 0. For (ii), note that since ttl : 
Ei^i — 7rL(Ei^i) is a diffeomorphism, and (x, ,^3) are coordinates on 7rL(Ei^i), 
it follows that (7r2(,T), 7r2(^3)) = {x,dx3S{x,yi,r]')) form coordinates on Ei^i. 
Since TEi = span{TEi 1, for both v — vr — 9-^ and v — vl — dy^, 
{x,yi,dx3S{x,yi, fjTjs)) are coordinates on Ei, and we can replace d^^S by 
?73 since dl.^^^S 7^ 0. Then, Eu is a hypersurface in Ei transverse to both 
dx^ and dy-^, and can thus be described as {r]2 — f{x,yi,r]3), yi — gixjij^)} 
with dx^g 7^ 0. 

Since {x,yi,r)') form microlocal coordinates on C, there exists a generat- 
ing function S{x, yi, 9'), so that y, 9') — S{x, yi, 9') —y' -9' parametrizes 
C near c°. We have 

C = {(a;, dxS- yi, de^S, -dy,S, 9') : x e M^ yi e M, ^' e \ o}, 

with 'KL{x,yi,9') = (x, dxS{x,yi,9')), ■KR{x,yi,9') = (y, -dy^S{x,yi,9'),9') 
and Ei,i = {92 - f ^ yi - g = 0}. Since 7rL(Ei,i) = {^1 = 6 = 0}, one 
has 9xiS'|{02-/=s/i-9=o} = and dx2S\{e2~f=yi-g=o} = 0. From the second 
relation, it follows that S'|{6i2_/=y^_g=o} is independent of X2, so that 5* is a 
function of just Xi, 0:3, 6'3; however, from the first relation one can express S 
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as 



S{x, y,, 9) = So{x^, 9^) + {y^ - gfS^{x, yi, 9') + 



f)S^{x,y,,9') 



Similarly, 7rfl(Sij) = {771 = 772 = 0} implies that dy^S\{e.,-f=Q=y^^g} = 0, 
but doesn't provide any more information about 5*, cf. Remark 1. On C, 
the ideal of Ei,i equals (6^2 - f{x,yi,93), yi - g{x,92.)) = {Sy^,92), so / = 
aSy^ + 1392. Hence /|si,i = 0; Sgg/ls^ ^ = (since / is homogeneous of degree 
1); / = c{yi — gY for some c e C°°; and de^g\Y,-^ ^ — 0, since g is homogeneous 
of degree 0. Thus, the canonical relation must have the form, 

C = \{xi,X2,X3,-2{yi- g)d:,^gSi + {yi- gfd.^^Si + {92- f)d^^S2- d^JS2, 

ivi - gfd^iSi - 2(1/1 - g)d^^gSi + {92 - 1)0^2^2 - d^JS2, 

{y, - gfd,,S, - 2(yi - g)d,,gSi + {92 - m,S2 - d,JS2 + ^,3^0; 



yi,y2,y3, '2{yi - g)Si + {yi - gfdy^Si + 



f)dy,S2-dyJS2, 



-O2, -93) : 2/2 = {yi - gfdg,S, + {92 - f)de,S2 + S2, 

y3 = (yi - gfOesSi - 2(yi - g)de,gSi + (^2 - f)d0,S2 - deJS2 + 9^3-^3 



From this, one sees that 



1 
1 
1 









-2d,,gSi 



9xiyifS2 
9x2yifS2 








dxiS2 
9x2^2 

8x382 











mod {yi-g,92-f) 



with det dTiL = [{-2dx^gSi-dx^yJS2)dx2S2-{-'2dx2gSi-dx2yJS2)dx^S2\{dl^Q^So 
+ {yi ~ 9)^1 + (6*2 — f)E2 for some Ei, E2. Since ttl has a cusp singularity 
at Ei^i, the determinant must vanish simply. We have that dl.^Q^SQ\Y,^ ^ 7^ 
by homogeneity and 8x2628 = c^a;2'S'2|si,i 7^ since dnL{TC) = {d^i = 0}. 
Since Ker dni = dy-^, the yi column in dni must be at Si 1. We have that 



= so (5'i+cS'2)|ei 



0, and 



{-2dx^gSi-d.^,yJS2)\E-,^, = dx-,g{Si+cS2 
{-2dx2g Si- dx2y J 82)1^:1,1 = dx29{Si+cS2)\T;r,i = 0, and since 5a;2S'2, 7^ 
we can conclude that det dTTilsi^i = ('S'i+cS'2) |si,i ■ Since So is nondegenerate, 
by another canonical transformation on the left we may assume that 5*0 = 
X393. From ker (iTTi — K-^, we have 5j/i(S'i+cS'2)|si,i — 0. Since d'KR{TC) — 
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{drji = 0} we obtain that 9j^i5'2|si,i = dg^yiSlsi^i = 0. One also has that 
= 0, since a calculation shows that this is dy_^rji{x, yi, 6'')|si,i, which 
vanishes due to the cusp structure of dirn and the fact that dT{R{TC) = 
{drji = 0}. Then one can write 5*1 = {yi — g^S^ + {62 — f)S% with 5*5 7^ 
0, leading to a weak normal form for a phase function parametrizing the 
flat two-sided cusp C, namely y, 62, O3) = (5*2 + {yi — g)'^SQ) {62 — f) + 
(a;3 - yz)0-i + {yi - gYS^ - y202, with d^^S2{x, y^ 6^') |si,i ^ 0, 5*5, 5*6 7^ 0. 

In studying the composition in the next section, /, g introduce only al- 
gebraic complications and do not affect the final result. Thus, for simplicity 
we take / = and g = xi and will use the phase function 

0(x,y,6'2,6'3) = (Sa - ^2 + - Z/l)^5'4) 6*2 + (X3 - ?/3)6'3 + (Xi - ?/i)'^S'3, 

d.,S2{x, yi, e') ^ 0, ^3, ^4 7^ on (9) 

with 5*5, Sq relabeled as S3, S4, resp. 



5 Composition 

Now consider the composition B*A for A G /™(X, Y; C), B G /™'(X, F; C). 
Conjugating by unitary FIOs associated with the canonical transformations 
used in the last section, we may assume that C is parametrized by the phase 
function (j) from (|9]). We have, for a G 5'™'+ 2, b G ^'"'"'"z, 

where the phase (f){z, y, rj') — (j){z, x, 6'), which we denote 0, is 

(2:3 - y3)Vi + (zi - yi)'^S3{z, yi,ri') + {S2iz, 77') - 1/2 + (zi - yifS^^z, yi, 77')) V2 
-(^3 - X3)6'3 - {zi - xiYS^iz, Xi, 0') - {S2{z, Xi,9') -X2 + {zi - XifSi{z, Xi, 9')) 62. 

We perform stationary phase in Z2,f]2 and z^^rj^: 

dz2$ = V2dz2S2 - 02dz2S2 + dz^Ssizi - yi)^ - d^^^Ssizi - xi)^ (10) 

+{zi - yifmdziSi - (zi - xif92dz2Si, 
dz-J) = r/3 - 6^3 + {zi - yiYdz^S'i - {zi - xi)'^dz3S-i + r]2dzsS2 - 02dz^S2 

+ {zi - yifr]2dz^Si - {zi - xif02dz^Si, (11) 
a^20 = {zi - yiYdr^^Sj, + 5*2 - 1/2 + {zi - y\)\S4 + r]2dr^^S^ + ?729,,2 5*2(12) 

a^30 = 2:3 - Z/S + {Zx - yxYdr^^S^ + ^2^ryiS2 + (^1 - ^/l )^r?2<9^3 'S'4- (13) 
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We solve for using ( 1T5]) and find 

= y3 - (zi - yi)^dr;^S3 - r]2dr,.,S2 - (zi - yiyr]2d^,^S4 := y3 + S^. 

From equation ( TT^ we can solve for Z2 since 82282 7^ 0, but only implicitly. 
Now rewrite ffTTl) as 



+ (zi - XiY92d;,^S4 - {zi - yifri2d:,^Si + 92d:,^S2 - r]2dziS2, 

and solve for 773, using the fact that 620^.^82 — r]2dz-jS2 vanishes at {xi = 
yi,r]' = 9'}, so that 

^2a.3^2-r]29.352 = (a:i-yi)e2a,23Xl^2 + (^2-r/2)(a.3^2+^2a,23e,52) + (^3-r/3)^29,%352. 

The other differences in the equations flTU]) . ([T5]) can be handled similarly, so 

^73 = ^'3 + (2:1 - yi)Mi{zi,Z2,Xi,yi,9',r]2) + (6*2 - V2)Ni{zi, Z2,Xi,yi,9',r]2) + 

(6*3 - T]s)P{zi, Z2, xi, yi, 9', r]2) 
for some functions Mi, A^^i and P. Solving for rj^ — 9^, we have 

V3-93 = [{xi - yi)Mi + {92 - r/2)iVi](/ + P)'^ = (xi - yi)M2 + (^2 - V2)N2. 

Also, r/2 - 6*2 = {xi - yi)M3 + (173 - 6'3)A^3, so that r/3 - 6*3 = {xi - yi)Mi and 
V2 — 92 = {xi — yi)N4,. Hence, the phase becomes 

= (2^3 - ^3)^3 + {X2 - y2)92 + ((2:1 - yi)^S4{Zi, y, 9') - {Zi - XifSi{zi, Xi, ?/2, 1/3, 9')) 92 

+ {zi - yiYSj, - {zi - X1YS3 + {S2{zi, y, 9') - S2{zi, xi, ^2, 1/3, ^')) ^2 + {xi - yi)SQ. 
We have 

{{zi-yifSi{zi,y,9') - {zi - xifS4{zi,xi,y2,y3,9'))92 

= [{zi - yiyS4{zi, yi, ■) - {zi - xiyS4{zi, yi, ■)]92 

+ [{zi - xifS4{zi,yi, ■) - {zi - xifS4{zi,xi, ■)]92 
= [{xi - yi){2zi - xi - yi)S4{zi,yi, ■) + (zi - - yi)dy^S4]92, 
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(zi-yiYSs- {zi- xiYSs = {zi-yiYSs{zi,yi,-) - {zi- xi)^Ss{zi,yi,-) 

+ {zi - Xi)^5'3(zi, yi, ■) - {Zi - Xi fS3{zi,Xi, ■) 

= (xi - yi){2zi - xi - yi){{zi - yif + {zi - Xiy)Ss{zi,yi, ■) 
+ {zi - xiY{xi - yi)dy^S3, 

and can write (zi — xi)^{xi — yi)dy^Si62 + {zi — XiY{xi — yi)dy^S2, + 
{S2{zi,y,9') - S2{zi,xi,y2,y3,0'))92 + {xi - yi)Se := (xi - yi)S7. From 
this, it follows that (j) = {x2 - y2)02 + (a^a - 1/3)^3 + (2^1 - Vi) [(2^1 - Xi - 
yi){{z\ - yif + {zi - Xif)S3 + {2zi -xi- yi)Si92 + 6*7(^1, y, Xi, 9)] . We make 
a change of variables similar to that in the model case in Sec. 12.41 

9i = {2zi -xi~ yi){{zi - yif + {zi - a;i)^)5'3 + {2zi - xi - yi)Si92 + 6*7, 
39 

= 2S,92 + {2z,-x^-yr)92d,,S, + 2S3{{z^-x^f + {z^-y{)'') 



^ = Ad,,S^92 + {2zi-xi-yi)92d\S^ + A{{z^-xif + {zi-yif)d,,S3 

+12(2zi -xi- yi)S3 + A{2zi - xi - yifd^.S^ + {2zi - xi - yi){{zi - yif 



+2(2^1 - xi - ^1)2^3 + {2z^ - xi - - y^f + {z^ - Xrf)d,,S3 + d,,Sj 



and 




implies that 



{xi - yi) 



^2S3 + Rd,,S3 



^Pi{xi,y,93) 



{xi - yif 



N{xi,y,93)+P,{xi,y,93) 



2 



2S4 + Rdzi 5*4 



2 
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for some N, Pi, and 

9i = P(i^i±^ + 2^2)^3 + RS,e2 + Sr := P2(xi, y, 9^). 

Since 5*3, 5*4 7^ and P = at the base point, one has 7^ 0. Following the 
same analysis as in the model case, one sees that 

(f>^{x2- y2)02 + {X3 - Z/3)^3 + (^^l - yi)^l + ( \ ' - ^1^2 - 7T^D<^1 + 
(73 Z ^3 ^3 

Next, we perform a stationary phase in t2,(J2- 

= -hai - (72, da^^= ^ + - yif - ^ - ^2 

(73 ^(73 (73 

and 

^ / A P2 A N , ,2 ^In 1 

= (x - y)e + (-i - - ti[-^ + — (xi - - -i] - -t?)ai 

C73 ^73 ^2, ^ 

With one more stationary phase in a\ and ^1, obtaining 
9010 = a;i - 7/1 + — , a^i0 = — - — - + — (xi - 7/1) - - 

t73 f73 ^3 ^3 Zfs ^3 Z 

and then homogenize by setting r — t\Q-i. The resulting phase is then 
^^{x2-y2 + j^{xi-yi)\e2 + (x^-y^ + ixi-yi) ' (14) 



^3 

+ {xl-ylf^^^ + {x^-y,)^{^fys, 
and the Schwartz kernel has the representation 

KpTix, y)^ j e^^c{x, y, 0', r) d02d0sdT, c e 3"'+"''. (15) 
For each of the variables s — xi, yi, y2, ys, let [[s]] denote 

r T N T T N 1 T 

^]]=^4^2 - ^Pi)+(a;i-2/i)V^3y, [N] = 5.3(^2-^Pi+(a;i-|/i)Vy+^^ 

{73 {73 Z (73 C73 Z Z 17 
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Then we have 

WF{F*F) C I (xi, X2, X3, ^ + 3(xi - yi) + P2 - ^Pi + ^(^)'^3 + (xi - yi)[[xi]ie'- 

L 6*3 6/3 6*3 

Z/i, 1/2, 2/3, ^ + 3(a;i - Vifr^ + - ^P, + ^(^)'^^3 + (xi - 

t/3 Zt)3 ^3 ^3 

6'2 + (xi - 2/1) [[2/2]], 6*3 + (xi -2/1) [[2/3 

= (xi - 2/1) (^^^ + 2(^) + - ^1)^ 

^ r 
= X2 -2/2 + -2/1);^ = 0, 

= 3:3 -2/3 - -yi)-Q^ + (^1 -2/1) [[^3]] = 0}. 

The phase function is degenerate because the equation dr4> = has a normal 
crossing; from {xi — yi = 0}, the contribution to WF{B*A) is contained in 

X,^ + P2-^Pl + ^(^f^3,^2,^3; 



x,^ + P2- ^Pi + O2, 03) --xe M^ (r, 0') e \ o|, 

[73 U3 Z t/3 J 

which is a cusp parametrization of A (as in the model case). On the other 
hand, if xi — ?/i 7^ then 



02 






= Pl- 


-(- 






y3 


= X3 - 


(Xi - 










= P2- 












= P2- 




V2 


= 02 + 


{xi - 



^2^ / ^^ 

) IT, 2/2 = a;2 + [Xi - 2/1)—, 
2 t/3 



2 



93 + riv ''"/''% (^i-yi)IN]. 

173 

h + rN ^""''^'^' + {xi-yi)[[yi]], 
[73 



Let C := be the image of ^^(xi, X2, X3, 6*3, r) = {x,^;y,ri) with ^2 = 
6*2,^3 = 6*3, and ^1, ?7i, 772, ^73, 2/2, 2/3 given above. As in the model case, ip 
satisfies Def. 12.21 it drops rank by one simply at Si = {xi — yi = = r}, 
and ker dvr = M ■ |: ^ TS, so that C is an open umbrella. 
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5.1 Generalized Fourier integral operators for A U C 

Using the phase 0, one can now define a class of generahzed FIOs, G, with 
the WF{G) contained in the union of the diagonal and the open umbrella. 
For a general a G S"^, define G by the right hand side of fllSI) . but with 
the amplitude a; then one has WF{G) C A U C. Furthermore, one can 
compute the orders of G on A \ C and C \ A, as in the model case at 
the end of Sec. 12. 4[ Due to the normal crossing of ^,-0, the critical set 
Crit-7 = {{x, y, r, 9') : d^-cp = 0, doicp = 0} decomposes as Grit/^ U Gritg, and 



where En 



D{x,e) 



and A,- are local coordinates on Grit- 



On Grit^, local coordinates are {x,t,6') and 



nf-T ^ Ql d(f> d(t> d(t> 

^y-^, ' ^ 862' de-i 



D{X,T, 6", 7/1,2/2, 2/3) 



02 
O3 



3r2 
2 61 



Pi 
O3 



hence G G /^(A \ C) and the principal symbol satisfies a ~ 5 2 , where 6 is 
the distance to A fl G. On the other hand, (x, yi, r, 6*3) are local coordinates 

on Gritc-, 



E, 



D{x,yi,T, 6*3, 



dr ' 96*2 ' de-i ' 



D{x,yi,T, 02,6s,y2,y3,O2 
2 



xi - yi 



+ 



9. 



G G I^{G \ A) and again a ~ 5~2. The equality of the orders on A and 
G away from A fl G is consistent with the composition result in |28l |6] for 
two-sided folds. This gives a precise description of the nonremovable artifact 
in the linearized seismic inversion problem for the single source geometry in 
the presence of cusp caustics. 
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